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UNIVERSAL NARROW-BAND EMI FILTER 

Abstract 

A new method of elimination of periodic! or quasi-periodic noise is presented. Present 
method employs, a« its main building block, a recently developed signal processing algo- 
rithm capable of extracting a specified compoiu^nt of a signal and tracking its variations 
over time. Performance of the present method is exemplified by application of the present 
EMI filter to elimination of power line noise potentially present on electrocardiogram and 
telephone cables. Superior performance of the method in terms of effective elimination of 
noise under frequency varying condition of pow<T line signal is observed. Present method 
offers a simple and robust structure which complies with practical constraints involved in 
the problem such as low computational resoura^s available and low sampling frequency. 
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Introduction 

Power line interference coupled to signal carrying cables is particularly troublesome in 
medical equipment such as electrocardiograms (ECGs). Cables carrying EGG signals from 
the examination room to the monitoring equipment are susceptible to electromagnetic 
interference (EMI) of power frequency (50 Hz or 60 Hz) by ubiquitous supply lines and 
plugs so much so that sometimes the EGG signal is totally masked by this type of noise. 
Filtering such EMI signal is a challenging problem given that the frequency of the time 
varying power line signal lies within the frequency range of the EGG signal. There are 
some other technical difficulties involved, the m(»st important of which is the low sampling 
frecjuency at which the EGG signals are taken and low computational resources available. 

The history of attempts to mitigate pow(T line EMI in EGG signals goes as far back as 
the EGG equipment itself. This problem was one of the first to attract the attention of 
developers of adaptive filtering theory [1]. Although classical adaptive filtering provides a 
partial solution to the problem, the problem is still considered open and research continues 
to find an ultimate solution [2, 3j. 

Pollution of EGG signals presents a general problem with medical equipment. This same 
problem may (K^cur in a variety of other scenarios. For example, telephone lines carrying 
voice or data are subject to indticed EMI from |)ower lines both in the form of differential 
mode and common mode interference. While elimination of connnon mode EMI is trivial, 
in practice some residual differential mode interfc^rence always persists to exist. Presence of 
such differtuitial mode EMI, frequency content of whic:h lies within the frequency spectrum 
of the signal of intt^rest, degrades the quality of the connriunication channel. The affected 
signal may be voice, or data in the case of telephone line-based data connnunications. 
The fact that characteriHtics of the interfering signal including its frequency may vary 
over time reuflers the noise suppression task difficult. Various methods of reduction of 
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power line iiiterfereiK:e have been proposed over the years each presenting strengths and 
weaknesses. No unique sohitiou to this seemingly simple problem lias been proposed so far 
[4]. Various signal processing schemes have been proposed in recent years for elimination 
of power line interference [3, 5], 

A recently developed signal processing algorithm, introduced in [6] was found promising 
in construction (if a universal EMI filter suitable for various applications in which the 
interference is a time-varying periodic (i.e. quasi-periodic) signal. It offers a robust 
structure and is shown to have a high degree of inununity with respect to external noise. 
This document presents structure and performance of such an adaptive EMI filter for the 
elimination of narrow-band interferences. Two examples of its application to ECG signals 
and signals carried by telephone lines are considered although the filter is general and 
may be applicable to other prol)lems of similar nature. 



Review of the Core Unit 

This section reviews the mathematical structun^ and i:>roperties of the core unit employed 
to construct the EMI filter of this invention. Let a(t) denote a signal comprising a number 
of individual sinusoidal components and noise, expressed by 

i/.(/)-^^fcsin<?^*.-fn(^) (1) 

k.=\ 

where 0a: = ^kt 4- 0^ is the total phase, and n{f) denotes the total noise imposed on the 
signal. The objective is to find a scheme for estimating a certain component of sudi input 
signal as fast and accurate as possible; a scheme which should not be sensitive to the noise 
and the potential time variations of the input signal. Simplicity of the structure, for the 
sake of practical feasibility, is desirable. 
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Let M be a manifold conbviiiing all pure siimsoidal signals defined as 

M = {y{t,e) = 9i Sm{02t 4- e [e.-.min.^i.max]} 

where 0 = |fli,<925^3]^ is ♦^he vector of parameters which belongs to the parameters space 

0 {[^1,^2; ^3)^ I ^ [^i.iiiiii' ^i.inax] } 

and T denotes matrix transposition. To extract a certain sinusoidal component, of the 
solution has to be an orthogonal projection of //(/) onto tlie manifold M, or equivalently 
has to be an optimum 0 whi(-li minimizes a distance function d between y(t,0) and u{t), 
i.e., 

^opi ^ argnunr/(/y(/,,ff), //.(O)- 
In least squares method, cl is the instantaneous distance function 

dHt0) = {u(t) -y{tj))'^c{t)\ 
The parameter vector 0 is estimated using the gradient descent method, 

where the positive diagonal matrix is the algorithm regulating constant. It controls the 
convergence rat<^ as well as the stability of the algorithm. 

Following the steps outUned above, a set of differential equations is obtained. The gov- 
erning set of equations of this algorithms can l)t* written as 

.4 = //-i<^sin<jl>, (2) 

uj — fi2eAcosd), (3) 

<p = M3ei4cos04-u;, (4) 

y{t) = Asm(l>, (5) 

e(0 = u{t)^y{t), (6) 
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Figure 1: Block diagram impleiiientation of the. core unit, 
where y(t) is the output. 

It has been shown that the ciynamical system n^presented by the above set of differential 
equations possesses a unique asymptotically stable periodic orbit which lies in a neighbor- 
hood of the orbit associated with the desired component of the function u{t). In terms of 
the engineering performance of the system, this indic^ates that the output of the system 
y{t) = Aiim<l> will approach a sinusoidal comi)onent of tlie input signal Ji{t). Moreover, 
time variations of parameters in u{t) are tolerat ed by the system. 

Figure 1 shows implementation of the algorithm in the form of c-oniposition of simple 
blocks suitable for schematic software development tools. Numerically, a possible way of 
writing the set of eq\iations governing th(^ presc^nt algorithm in discrete form, which can 
be readily used in any progrannning language, is 

A[n -1-1] = A[rt] -4- Tsfhe\n] sin 0[7/.], 
Lj[n H- 1] = uj[n] + r^//2^Mi4[n] cos (t>[n , 
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(l)[n + 1] = (l)[n\ -h T,uj[n] + T^mi:ie[n]A[n] cos (^[nj, 
y[n] = A[n]siu0[n], 
e[n] = '«[n] — i/fn] 

where a first order approximation for derivatives is assumed, is the sampUng time and 
n is the index of iteration. 

In the simulations presented in this document, Matlab SiauUink^*^'^ computational soft- 
ware is used as the main computational tool. Figure 2 shows an example of performance 
of the core algorithm in which the frequency of the input signal undergoes a step change 
of 10%. It is observed that the variations are effectively tracked with a transient of few 
cycles. Values of the parameters are chosen to be //i = 100, //,2 = 10000, //;^ = 0.02. for 
this simulation. 

The dynamics of the core algorithm presents a notch filter in the sense that it extracts 
(i.e. lets pass through) one specific sinusoidal component and rejects all other components 
including noise. It is adaptive in the sense that the notch filter accommodates variations of 
the characteristics of the desired output over time. The center frequency of such adaptive 
notch filter is specified by the initial condition of frequency uj. In Figure 1 this initial 
value is explicitly shown for easy visualization. 
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ire. 2: Response of the algorithm to a step change in the frequency of tlie input signal. 
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Structure of the Present EMI Filter 

One single unit of the core algorithm can be t^niployed to extract the quasi-periodic 
interference mixed with the signal. This unit can effectively follow time variations in 
amplitude, phase and frequency of the interfering signal. Once it is extracted, it can be 
subtracted from the input signal to yield a clean signal. 

In order to improve performance of the unit, the use of a band pass filter (BPF) to filter 
out non-interference signal components is i>roi)Osed in Figure 3. This Ijand pass filter 
does not need to be sophisticated and can be as simple as a second order filter. Its role 
is to improve the signal to noise ratio (signal here meaning the interference and noise 
meaning all other components) at the input of the core unit. Whatever is not removed 
by this BPF will be cflFectively removed by the core unit so as to produce a single pure 
sinusoid which is the interference. This interference in then to be subtracted from the 
input to provide the clean signal. BPF characterized by its transfer function H{f) causes 
an attenuation \H{f)\ and phase delay ///(/) which are functions of frequency. Since 
the core unit generates the value of the frequency in real-time, the attenuation and phiise 
delays are known and can l)e restored as depicted in Figure 3. In general, the filter does 
not have to l^e band pass and a high pass filter may l)e sufficient. As a concrete example, 
the band- pass filter employed as the pre-filtering tool in Figures 3 and 4 is chosen to be 
a second order filter characterized i)y the following transfer function: 

_ l()0<s 

Gain and phase characteristics of this filter are shown in Figure 5 in which 'jJo is taken to 
be 100 Hz for the ease of visualization. 

Where the interfering signal is severely distorted, the harmonics may also be present. In 
such cases, desired signal is not polluted by a sinusoid, but by a number of sine waves. A 
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Figiiro 3: Block diagram of the present EMI filter. 
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Figure 4: A configuration to eliminate a severely distorted quasi-periodic EMI. 
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Figure 5: Gain and pluuso diagrams of the pro-filter. 

more general coiifiguration in Figure 4 may then be employed to eliminate the fundamental 
and the harmonics of the EMI. 

Application of the Present Filter to ECG Signals 

ECG signal is basically an index of the functionality of heart. The physician can detect 
arrhythmia by studying abnormahties in tlie ECG signal. Since delicacies present in the 
ECG signal convey important information, it is important to have the signal as clean as 
possible. Figure 6 shows a clean ECG signal recorder! at Beth Israel Hospital (BIT) in 
Boston and made available by Massachuset Institute of Technology MIT-BIH [7]. The 
recording was done using battery operated EGG equipment to minimize power line EMI 
although some such EMI still exists which is mostly coupled at the time of recording 
the signals on the tape. The frequency spectrum of this signal spans from near DC 
frequencies to about 100 Hz. The sampling frecjuency in most ECG devices are 240 Hz 
or 360 Hz. In this case^ the equii)ment was operated at the sampling rate of 3G0 Hz. 
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Figure 6: Recorded ECG signal and its frequency spectrum. 

Therefore, the spectrum can theoretically include frequencies from zero to 180 Hz. The 
data available from MIT-BIH contains a liigh DC offset wliich is eliminated in sinuilations 
of this document. 

ECG signals can be easily polluted by power line noise of relatively large amplitude. Were 
the frequency of power line interference acxuratt^ly at 50 Hz or 60 Hz, a sharp notch filter 
would be able to separate and eliminate the noise [8]. The major difficulty is that the 
frequency can vary about fractions of a Hertz, or even a few Hertz in some countries. The 
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sharper the notch filter is designed, the more inoperative, or rather destructive, it becomes 
if any change in the frequency of the power Hne occurs. Of course, turning the notch 
filter into a band stop filter by widening its rej(^ction band, and thereby acconnnodating 
frequency variations, does not offer any better solution since it will undesirably distort 
the ECG signal itself. The freciuoncy of the power grid is usually taken as being constant 
when conventional EMI filters for ECGs are designed. In such arrangements, the system 
is very fragile with respect to power frequency variations and can become completely 
malfunctioned. SucJi adaptive or non-adaptive filters, those discussed in [9] for instance, 
greatly suffer from this shortcoming. 

One of the possible alternatives to take frequency variations into account is the use of 
external reference power line signal [10]. This t<T.hnic|ue, available by the use of adaptive 
filters only, is reported practically difficult or rather impossible [9], For this reason, other 
methods, usually very complex and inflexible, are constantly being proposed [2, 3]. 

The ideal EMI filter for ECG is the one which act« as a sharp notch filter to eliminate only 
the undesirable power line interference while automatically adapting itself to variations 
in the frequency of the noise. Of course, this adaptation must be done very quickly 
so as to keep the signal clean all the time. It is supposed to he able to work in low 
information background, namely that dictated by low sampling frequency, and must be 
robust with respect to variations in its internal as well as external conditions. An example 
of internal condition is its settings. External conditions can range from the temperature 
of the environment in which the equipment is supposed to function to the superimposed 
noise/dis tort ion on the interfering power signal. 
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Figure 7: A power line interference of 1 mv level at 60 Hz is added to the ecg signal to 
provide a i)ollnted input signal for the present filter. 
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Figure 8: Performance of the filter in eliminating a power line interference of 1 mv level 
at GO Hz. 
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Figure 9: Porforinance of the filter in eliininatiiig power line interfereuc^e of time varying 

level (oscillating between 0 to 1 niv) at 60 Hz. 
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Figure 10: Performance of the filter in eliminating power line interference of 1 mv level 
at 55 Hz. 
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Figure 11: Performance of the filter in eliminating power line interference of time varying 
level at C5 Hz. 
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Performance of Present EMI Filter in the Case of ECG Signals 

This section presents a number of simulations to demonstrate performance of an EMI 
filter constructed as in Figure 3. The filter at the input of the core unit attenuates the 
ECG signal by a factor of 10. Figures 7 and 8 show performance of the filter in eliminating 
a power line interference of 1 mv constant level whose frequency is fixed at 60 Hz. This 
is an elementary experiment since a simple notch filter will easily eliminate such a fixed 
frequency noise, hi the sinmlation, only a 60 Hz component is added to the clean ECG 
signal obtained from MIT-BIH [7] whose presence is clear in the frequency spectrum of 
the input (Figure 8). 

The values of the parameters /ii, /i-j, and /i3 determine the convergence speed versus error 
compromise. Generally, the higher the values are chosen, the faster the algorithm tracks 
variations at the expense of larger steady state <^rror. Therefore, it is important to define 
desirable balance of speed/error. For the simulations in this document, a moderate choice 
of fi] = 200, fji2 = 20000, and /i3 = 0.01 results in an EMI reduction of a factor of about 
20 while keeping the transient time short enough. All initial conditions of integrators are 
zero except for that of frequency which is taken to be 60. 

To demonstrate the ability of the filter in adaptively tracking the variations of the level 
of noise, the level of the EMI is made to be changing with time in Figure 9. Again, error 
is confined to within about ^ of the maximum noise in the input. 

Figure 10 shows performance of the filter in tracking the variations in frequency of the 
power line noise. The filter is adjusted -by virtue of its initial conditions- to extract a 
power line noise of frequency 60 Hz. The EMI in the input, however, is oscillating at 
55Hz. Effective tnicking of unknown input frequency is observed. 
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Time Domain Representation (Normalized) 




Figure 12: Tiiiio-doniain roprosontation of the porfonnance of the present EMI eHminator 
in suppression of a fixed frequency interference from an arbitrary input signal. 

Finally, Figure 11 shows performance of the filter when all characteristics of the EMI 
are changing with time. As before, the filter is adjusted to extract a power fine noise of 
60 Hz. However, the incoming EMI has a frequency of 05 Hz. The level of the EMI is 
also changing with time. Again, eflFective elimination of superimposed EMI is observed. 



Performance of Present EMI Eliminator in the Case of Telephone 
Line Signals 

Figure 12 shows performance of the EMI filter in suppression of an interference of fixed 
frequency of 60 Hz from an arbitrary input signal. The original input signal is arbitrarily 
taken as a recorded bird chirp sampled at = 8192 Hz. It goes without saying that 
the notch filter targets the pre-specified sinusoidal interference present in its input, hence 
frequency (.omposition of the incoming signal is irrelevant as far as EMI filter is concerned. 
The level of interference added to the original signal is taken to be equal to that of the 
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Figure 13: Error incurred in suppression of a fixed frequency interference from an arbitrary 
input signal. 
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Figure 14: Froquency-doniain representation of the ])erf(;rniance of the preyeiit EMI eHm- 
inator in suppression of a fixed frequency intorf<Tenre from an arbitrary input signal. 
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Figure 15: Magnified spectrum of tlie original aud retrieved signals around the frequency 
of the interfering sinusoi<l (60 Hz). 

original signal so that signal to noise ratio (SNR) is 0 dB in the following simulations. 
Figure 13 shows the incurred error in removing the EMI. It is basically the difference 
between the refined signal and the original signal not yet polluted by the interference. 
Frequency spectrunis of the original, polluted and retrieval signals are shown in Figure 
14. Figure 15 shows the magnified portion of the frequency spectrum of the original and 
the retrieved signals about the frequency of the interference. It is observed that the notch 
filter effectively removes the interference and retains the original signal almost untouched. 

To show the adaptive nature of the present interference eliminator with respect to vari- 
ations in frequency of the interfering signal, performance of the present method when a 
step change in the frequency of interference occurs is shown in Figure IG. As noted be- 
fore, the rate of convergence in tracking time-variations such as that shown in Figure 16 
is totally controllable by means of adjustment of parameters- The tracking capability of 
the present method is advantageous in devising a universal power supply noise eliminator 
nidependent of the nominal frequency of the grid. Regardless of initial frequency setting 
(whether 50 Hz or 60 Hz), the filter finds the instantaneous frequency of the power line 
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Frequency Change 




Time (8) 

Figure 16: Response of the present method to a frequency step change in the interfering 
sinusoid. Top figure shows the estimated frequency and the bottom figure shows the 
incurred error in EMI suppression. 

Frequttncy Offset 




Time (s) 

Figure 17: Response of the present method to a large frequency offset in the interfering 
sinusoid. Top figure shows the estimated frequency and the bottom figure shows the 
retrieved signal. 
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Figure 18: Performance of the present method in eliminating a sinusoidal interference 
of time-varying amplitude. The top figure is the interference, the middle figure is the 
estimated value of the amplitude of the interference and the bottom figure is the incurred 
error in EMI suppression. 

interference. Figure 17 shows an example in whicli the EMI filter expects a 60 Hz power 
supply interference while the interference happens to be of 50 Hz frequency. 

To demonstrate adaptive nature of the present method in tracking variations in the level 
of the interference, performance of tlie present method in eliminating an interference of 
time-varying amplitude is shown in Figure 18. 

As noted before, power supply interference may be distorte{l by the presence of harmonics. 
Figure 19 sliows an example of a highly distorted interference in which harmonics of third 
and seventh order are present. In such cases, a multiplicity of EMI filters, connected in 
parallel combination as suggested by Figure 4, may be used for EMI mitigation. Figure 20 
shows the frequency domain representation of the performajice of an EMI filter consi.sting 
of three units connected in parallel. The fundamental frequency of the interference is 
made to slightly vary over time. An EMI reduction of about 40 dB is observed. Again, 
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Figure 19: Interference is taken to he highly distorted by the presence of harmonics. This 
figure shows one cycle of the interfering signal. 
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Figure 20: Frequency domain representation of the performance of the present method in 
elimination of interference of liighly distorted shape. An EMI reduction of about 40 dB 
is observed. 
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the level of desired EMI mitigation is controllable by the adjustment of parameters and 
at the expense of speed. Considering that the power supply noise is usually a slowly 
time-varying signal, one can sacrifice speed for better interference elimination. 
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